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Abstract 



We have investigated edge modes of different multipolarity sustained by 
quantum dots submitted to external magnetic fields. We present a microscopic 
description based on a variational solution of the equation of motion for any 
axially symmetric confining potential and multipole mode. Numerical results 
for dots with different number of electrons whose ground-state is described 
within a local Current Density Functional Theory are discussed. Two sum 
rules, which are exact within this theory, are derived. In the limit of a large 
neutral dot at B=0, we have shown that the classical hydrodynamic dispersion 



law for edge waves u){q) ~ \/q In (qo/q) holds when quantum and finite size 

effects are taken into account. 
PACS 73.20.Dx, 72.15.Rn, 78.20.Bh 
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I. INTRODUCTION 



Collective excitations induced in finite fermion systems by external probes have been 
extensively studied in the last years. Particular effort has been devoted to the understanding 
of the giant dipole resonance in nuclei and of the plasmon mode in metallic clusters 
0-^. Recently, a strong collective state has also been observed in quantum dot structures 
I^J^. These collective modes have all in common the feature of being L=l, S=0 excited 
states induced by a dipolar external radiation which is the dominant component of the 
electromagnetic field when its wavelength is much larger than the size of the system. Dipole 
modes correspond to oscillations of protons against neutrons in nuclei, and of electrons 
against the positive background in metal clusters and dot structures. Whereas in nuclei the 
restoring force of the collective motion is the symmetry potential which acts differently on 
protons than on neutrons, in metal clusters and dots it arises from the Coulomb interaction 
between ions and electrons. 

Collective states with multipolarity L>1 were predicted to exist and have been observed 
in nuclei. At present, we have a complete systematics of quadrupole (L=2) and octupole 
(L=3) collective excitations in nuclei 0,0,3 • Multipole collective states in '^He droplets have 
also been studied P-[Tl[|, but have not been experimentally detected so far. In metal clusters, 
the predicted multipole states |T^,|T3[ have not been observed either. The basic reason is the 
experimental difficulties arising from the fact that free clusters are produced and analyzed 
on-fly, with the added difficulty in the case of ^He drops of being electrically neutral systems. 

Far-infrared absorption spectroscopy experiments on large radius quantum dots submit- 
ted to a static external magnetic field B have likely evidenced quadrupole excitations 0, and 
an anticrossing between the L=l and L=2 resonances, each of them splitted in two branches, 
one with negative and another with positive B-dispersion. We recall that for these systems, 
L has to be understood as the angular momentum about an axis perpendicular to the dot 
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plane. 

The collective spectrum of quantum dots has been addressed in different ways. An 
explanation based on classical arguments can be found in p|. In Refs. [^^-|T6|, use has been 
made of an edge-magnetoplasmon model, whereas in Refs. |T7| , p^ a Hartree-RPA method 
has been employed, and a Hartree-Fock-RPA one in Ref. [M. Multipole modes have been 
considered in Refs. and JlSl . All these approaches neglect the electron correlation energy, 
and therefore, no microscopic calculation of L>1 modes exists which incorporates exchange 
correlation energy as well as quantum effects arising from finite size and shell structure of 
dots in a magnetic field. A workable, yet reliable formalism that takes into account all 
these effects is called for to describe these excitations. In this paper we present one of such 
possible formalisms. It is based on the equation of motion method within the framework 
of a local Current Density Functional Theory (CDFT) [^. This functional theory is well 
suited to study electronic systems in presence of a magnetic field, and has been successfully 



employed to study ground state (g.s.) properties of quantum dots El 



This paper is organized as follows. We describe the equation of motion method in 
Section II, as well as the basics of the strength function and its moments (sum rules) needed 
to interpret the experimental and theoretical results. In Section III we present a rigorous 
solution of the dipole mode in the case of a parabolic lateral confining potential. The 
exactness of the dipole solution in this case, irrespective of the value of the magnetic field, 
has been discussed in [^,^]. Here, we have obtained this results in a way that yields an 
explicit expresion not only for the spectrum, but also for the eigenstates. In Section IV 
we present a variational approach to the description of multipole excitations. We show in 
Section V that for neutral large dots at B=0, these excitacions build an edge wave with 



dispersion relation of the type uj{q) ~ y q\n{qo/q). Detailed numerical results are presented 
in Section VI for dots with different number of electrons, and the concluding remarks are 
presented in Section VII. Finally, an appendix contains some technical details. 
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II. THE EQUATION OF MOTION METHOD 



We consider N electrons moving in the z=0 plane where they are confined by the dot 
potential V+(r) with r = a/x^ + y'^. On this system it may act a constant magnetic field in 
the z-direction described by the vector potential A = i(— y,x,0)B, and we suppose that the 
system can be described by the N-electron Pauli's Hamiltonian: 



N 



i=l 



e 



c 



2 I AN 



i<j 



2^ 

where m* is the electron effective mass which together with a dielectric constant e and 
gyromagnetic factor g* are characteristics of the semiconductor (for example, g* = —0.44, 
e=12.4 and m*=0.067 m^ in GaAs), fi*^ is the effective Bohr magneton fi*^ = he/{2m*c), a 
is the Pauli matrix vector, and Vc is the electron-electron (e-e) interaction 

t \ 1 1 I j \ 

Eq. can be rewritten as 

^ f »2 1 1 1 1 ^ 

H = E\^^ + ^^ch. + -m*u:lr1 + -g^^^^Bo^^ + \4(r,)| + ^ K(| - |) , (3) 

where a;c=eB/m*c is the cyclotron frequency and I2- is the angular momentum operator 
about the z-axis: 

= -'4 ■ 

Given the exact g.s. |0) of the N-electron system, it is possible to obtain the exact 
spectrum corresponding to a broad class of collective vibrations if one is able to find an 
operator 0+ such that the following equation of motion is fulfilled: 

\E, 0+] = TiM^ . (5) 

The state O"'"|0) has an excitation energy hoj, and the g.s. fulfills O|0)=0. As the excited 
states have a well defined angular momentum, so the operators O"^ must have. Consequently, 
one has to solve Eq. (|^) for each L- value. 



When a magnetic field acts perpendicularly on the dot, it causes a splitting of the excited 
B=0 states into two branches of energy hu)±L, each carrying an angular momemtum ±hL 
over that of the g.s. It implies that besides Eq. (||), the physically acceptable O^^ operators 
have to fulfill 

[L,,0+J = ±^LO+i (6) 

with L=l,2,3,... and = J2iLi hi- If ^Lq is the angular momentum of |0), the states 
I ±L) = I 0) have an angular momentum h(LQ± L). 

To easy the formulae that otherwise would be very cumbersome, from now on we shall be 
using effective atomic units, defined by ^ = e^/e = m* = 1. In this units system, the length 
unit is the Bohr radius ao times e/m*, and the energy unit is the Hartree times m*/e^, which 
we call respectively, ag and EJ^. For GaAs we have ~ 97.94 A and E^ ~11.86 meV. 

It is obvious that to find these operators is in general as difficult a task as to solve the 
Schrodinger equation corresponding to the Hamiltonian Eq. (|I]) for the vibrational states, 
and one is led to solve Eq. in an approximate way, one of this being, for example, 
the Random Phase Approximation. Another fruitful approximation, originally proposed by 
Feynman to describe density excitations of superfluid ^He , consists of making an ansatz 
on 0+1 0). Acting upon |0) with Eq. (|^) and projecting onto 0+ |0) one gets 



(0 


\[0,[H,0^]]\ 


0) 


(0 


[0,0+] 


0) 



(7) 



Eqs. @ and (|^) are completely equivalent if |0) is the exact ground state and 0+ is the 
sought operator. The advantage of Eq. is that we may look for approximate solutions 
of variational type, guessing 0+ and obtaining |0) within a workable, yet accurate scheme, 
such as the Local Density Approximation (LDA) at B=0, or CDFT at 0. To find 
these approximate solutions is the subject of Section IV. We shall also see that, remarkably, 
Eqs. (|1]) and (|^ have exact solutions in the dipole case when the confining potential has a 
parabolic form. 

The excitation spectrum of the system is usually probed by different external fields, or 
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given the well defined angular momentum of the excited states, by a selected L-polar compo- 
nent of the field. For an excitation operator F representing it, a useful, often experimentally 
accesible quantity, is the so-called strength function: 

S{E) =J2\{n\F\0)\'5{E- K) , (8) 

n 

where E„ and |n) are, respectively, the excitation energy and the excited state, and the sum 
or integral in the case of continuum spectrum extends over all excited states of the system. 
Of special interest are some energy moments of the strength function 

m,= [ E'S{E) = E^nl (n|F|0) f (9) 

which we shall call sum rules (SR). They are the mi and moments, which can be also 
written as 

mi = ^(0| [F+,[H,F]]\0) 

(10) 

ms=^-{0\[[H, [H,F+]],[H,F]]\0) . 

These SR have been extensively studied in the literature For the present purposes it 

is enough to recall that, if only one excited state is contributing to S(E), E3 = (ma/mi)^/^ 
coincides with the corresponding excitation energy. In a more physical situation, whenever 
the strength is concentrated in a narrow, "high" energy region, E3 is a fair approximation to 
the resonant energy. We call E3 the scaling energy because m3 can also be obtained scaling 
the |0) wave function as 

I V) = e'^""'^^ I 0) , (11) 
and then carrying out the following derivative: 

1 
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The ms sum rule can be obtained either from Eq. (|TOD or from Eq. (p!2D. The latter 
allows one to identify as the restoring force that causes the collective vibration, and 
consequently, to identify mi with the collective mass parameter. 

We want to stress that not all possible approximation schemes fulfill the mi and ms sum 
rules in the sense that a direct evaluation of Eq. (^) yields the same result as Eqs. ([I^). A 
test on the consistency of the approximations made to get the g.s. and excitation spectrum, 
is the fulfillment of these SR, especially of mi which is very model independent. 

III. EXACT SOLUTION OF THE DIPOLE MODE FOR THE PARABOLIC 

CONFINING POTENTIAL 

It has been frequently argued in the literature that V+(r)= |c<jQr^ can be a good ap- 
proximation to the confining potential. That is the case, for example, when the number of 
electrons in the dot is small as compared with the number of positive ions N+ that pro- 
duce the confining potential |^. If that approximation holds, it is easy to check that the 
solutions of Eq. (|^) for the Hamiltonian Eq. (j^) and L=l are 

(13) 

where 

N N 
i=l i=l 

(14) 

N N 
i=l i=l 

and 

u = . (15) 
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The corresponding frequencies are 

^± = ^ ± Y , (16) 

and it is easy to verify that 

[L„0+] = ±0+. (17) 

The states |0) are normahzed to unity and orthogonal, carrying an angular momentum 
(Lo±l). This exact result stems from the translational invariance of the electron-electron 



interaction for which 



Y.i<j Kd ri — fj \) , P =0, and consequently it is independent of the 



e-e interaction provided it is local [23 



It is a simple matter to check that the states | ±1) exhaust the mi and m^, SR for 
the dipole operator ^=YliLiXi, and that the dipole strength is equally distributed between 
them. One gets 

mi(D) = ^N (18) 

m,{D) = \ N{u' + (19) 

I (0 I D I +1) p=| (0 I I -1) |2= ~ , (20) 

and for example, 

CU+ I (0 I D I +1) |2 + cu. I (0 I D I -1) |2= ^ = mi(D) (21) 

It is interesting to notice that in the B=0 limit, = uj^= E^iT)) = ujq. This result is 
independent of the number of electrons in the dot, in agreement with the generalized Kohn 
theorem |2^ . The parabolicity of the potential is expected to break down when the number 
of the electrons in the dot increases and the electronic density extends up to the edge of the 
dot. Departure of V+(r) from the parabolic law originates an N-dependence in uj± and in 
E3(D). 
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The exact solution to Eq. can be used to obtain the static dipole polarizabihty of 
the dot a(D), which is twice the m_i(D) sum rule |1|J 



m_i(D) = — I (0 I D I +1) p +— I (0 I D I -1) \'= . (22) 

Consequently, in the parabolic potential approximation the dipole polarizabihty is indepen- 
dent of the magnetic field. We expect a{D) to be B-dependent in the case of a more general 
confining potential. Using that 

= ^ (23) 



we get 



a{D) = R'^ . (24) 



The R-dependence could have been anticipated from a dimensional analysis. 

Let us finally discuss the nature of the dipole modes with regard to their geometrical 
shape. If we consider the transition density associated with the excited states | ibl) 

p,,(f) = (0 I pO+ I 0) = (0 I [/5,0+] I 0) , (25) 

where p is the N-electron density operator 

TV 

p = J2s{f-f^, (26) 

1=1 

it is easy to obtain 

pUf)ocp'ir)e^'' , (27) 

where p(r) is the g.s. electronic density and the prime denotes the r-derivative. This 
transition density is characteristic of an edge excitation. 

IV. A VARIATIONAL SOLUTION FOR GENERAL AXIALLY SYMMETRIC 

POTENTIALS 



When the confining potential is no longer parabolic, irrespective of the value of the 
magnetic field and of the multipolarity of the mode, the equation of motion (^ cannot be 
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solved exactly, and we have resorted to an approximate method based on Eq. (|^). In the 
dipole case, the natural guess for the operators is: 

OX = a+{Q - ib+P) 

(28) 

0+ = a-{Q+ - ib-P+) , 

where b-t are variational parameters and a± have to be determined from the normalization 
condition. This process can indeed be carried out, and one gets the dipole spectrum of the 
system. 

At this point, we have considered it more convenient to generalize the Q and P operators 
in such a way that the calculation can be done for any L-value. To this end, we have first 
taken the following Q^, 

N N 

QL = T.^fe^''' =E?1 • (29) 

i=l 1=1 

This choice is inspired in that (kr)^ e*^^ is the small k-expansion of the function JL(kr) e'^^, 
which is the restriction to the z=0 plane of the general solution of the Laplace equation in 
cylindrical coordinates. Jl is the L-Bessel function of the first kind | p7|| . 

Next, we have taken as partner of in the O"*" expression the following operator: 

TV 

Rl = Y.^l-iP, (30) 

i=l 

where the i-particle index is implicit in the particle and momentum coordinates. This choice 
is again guided by the exact dipole case, since the conmutator [H, Ql] yields a combination 
of Ql and R^, which reduces to the one we have found in Section III when L=l. 
We have thus considered as O^^ operators the following combinations: 

OXl = a+iQi - ib+LRi) 

(31) 

OtL = a4Qt-^b^LRt) 
10 



Eqs. (|31|) reduce to Eqs. ( ]28|) for L=l, and the proposed operators fulfill Eq. (H). These 
modes are edge modes whose transition densities read 

PtHrl = (0 I [p,0+J I 0) (X r^-V'(r)e±^^^ . (32) 

To minimize Eq. (0), it is necessary to have a description of the g.s. |0). We have taken 
for it the Kohn-Sham sate built from single particle (s.p.) wave functions obtained in the 



framework of the CDFT of |2l[. We refer the reader to that work for the details. For the 
present purposes, it is enough to recall that the s.p. wave functions 4>nia{r,9) are separable 
in r and 6 variables as 

<PnUr,e) = Unia{r)e-''' , (33) 

with I = 0, ±1,±2, ... being the orbital angular momentum of the s.p. state. Upon mini- 
mization one gets: 



46±(mi ±6±L2r) 



(35) 



where 



ma , L 



mi 4 

r 1 



+ T^c (36) 



Ui = - 

mi mi 



= —(L-l) I dfr^'Sir) (37) 
mi J 

mi = L^ j dfr^^-^p{r) (38) 

"^3 = "^3(^) + "^3(ee) + m3(+e) (39) 
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with 



m3(T) = L^{L - 1) y dfr^^-^ [Lr + 2{L - 2)A] 



(40) 



2L 



2^L\ JO 



Er 



1 

r Jo 
r 



oL+i ri 



(2L + 1)!! 



,2L-3 



- 27rL^ / drp'{r) {r 

Jo JO 



2L 



El ( - ) rfr' - ^r2V'(r) 



m3(ee) = 471^^^^^^^ p'{r)dr \ [ \2{L + l)p'(r')r'2^ + p"(r')r'2^+'l ( ^ ) rfr' 



(41) 



m3(+e) = -TiL' I drV+{r) (2L - l)r"^-V(r) + r^^'V" 



Besides the particle and kinetic energy densities 



(42) 



p(r) 



nl(7 



(43) 



r(r) = (0 I E V 5(r - rl) vl 0) = E«i + -^1.) 



nla 



we have introduced in these equations the densities 



(44) 



li.'r) = - Y.^^lea{r) = rjp{r 



(45) 



A(r)=E 



as well as the E^-f unction: 



TT 



(46) 



(47) 



where F is the hypergeometrical function P^. The density A(r) represents a sort of centrifu- 
gal kinetic energy density. It can be shown that in the zeroth-order Thomas-Fermi (TF) 
approximation A(r)=|r(r), similar to the A (r) = |r(r) relation holding in three dimensions 
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28| . The function jp(r) is the paramagnetic current density It is worth to see that at 



B=0, 7(r) vanishes due to time reversal invariance. Its contribution at high magnetic fields 
is crucial to have well behaved B-dependent collective energies. Eqs. ( p^ - |4^ ) constitute the 
main outcome of the present Section and will be used in Section VI to obtain numerical 
results within CDFT. 

The goal of describing the multipole modes in a situation as general as possible makes 
these expressions to look very cumbersome at first sight. However, it can be noticed that: 

i) For L=l and the parabolic V+, one recovers the exact solution. 

ii) For any L and zero magnetic field, b± = 1/0)3 and O^i^ |0) merge into a degenerate 
state of energy 



u^±L = J— . (48) 
V mi 

This is precisely the B=0 scaling energy E^^Ql), since mi is actually the general mi^Qi) 
sum rule, and as it is outlined in the Appendix, ms reduces to the zero magnetic field m3(Qx,) 
sum rule making 7(r)=0. This result is the variational analog of the case we have discussed 
for the dipole mode in a parabolic V+. 

iii) In the L=l case, irrespective of B all terms in but rh^l+e) vanish, and for any 
axially symmetric ionic potential V+(r) we find 



^±1 = y ^ + ^ / drAV+ir) p{r) ± , (49) 

which reduces again to the exact case if V4. = \uJqt'^. 

iv) For any B, the states | ±L) and their energies satisfy the mi sum rule, i.e.: 

mi(Qi) =L'J dfr'^^-'pir) = 00+ \ {0 \ [Ql,0+J | 0) ^ + | (0 | [Ql, J | 0) 1^ . (50) 

It implies that, although O^^ |0) are only approximate solutions to the L-polar excitation 
spectrum, within our method there are no other excited collective states corresponding to 
that multipolarity built as coherent superposition of one electron-hole states. The sum rule 

3 

msiQL) = fns{T) + -L^cj^mi + StUcL^T + m^{ee) + m3(+e) (51) 
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is also fulfilled. Moreover, the square transition amplitudes to the | ±L) states are equal, 



{0\[Ql,OX^]\0) \' = \ {0\[Ql,O+^]\0) 



(52) 



The fulfillment of mi and constitutes a rather stringent test on our variational method. 

v) It is worth to notice that in all cases, there is no explicit contribution to these formulas 
arising from the exchange correlation energy terms in the CDFT Hamiltonian. This is 
because Ql is a solution of the two-dimensional Laplace equation. A similar result holds in 
three dimensions within LDA if one takes for a solution of the corresponding Laplace 
equation p9 |. 



V. EDGE WAVES IN NEUTRAL LARGE DOTS AT B=0 

It is possible to use the scaling energy E3, which is a good approximation to the collective 
excitation energy at B=0, to obtain the dispersion relation of the edge waves in the case of 
a neutral large dot. 

Under these conditions, the g.s. electronic density is constant everywhere apart from a 
narrow region along the border of the disk. Let po=N/(vrR^) be that density, and let a be 
the width of the edge region. Using the zeroth-order TF approximation to the kinetic and 
centrifugal kinetic energy densities 



one gets from Eqs. (p8| , ^0| -| 



r = np 



(53) 



mi = tipqLR 



2L 



(54) 



2 „2 d2L-2 



ms{T)=2L'{L-l)TT'pf,Ii 



(55) 
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m3(ee) + m^{+e) = At.L'R^''-^ pl ^'^^^r }}'' Fl{z) 



with z=l — (9(a/R) at least, and 

Fl{z) -= 



1f(-,L + -;L + 1-z' 
2 V2' 2' ' 



(56) 



(57) 



which diverges at z=l. Physically, this divergency is associated with the known divergency 
of the electric field at the edge of the disk in the case of a constant electron density ||3^, see 
below. 
Thus, 

L(L-l) , oo(2L-l)!! 



Ei = 27rpc 



+ 4^ 



:Fr(z) 



i?2 R2^{L-1)V 
Fi(z) can be written as a function of Fo(z), which is the elliptic function K(z) |]27[| . 



(58) 



(2L- 1)!! 
2^L! 



1 



Using it we get 



L(L-l) , L 



\2m-l 



1 ^1 

lfjz) + 1 - y —— 

2 ' i'i2m-l 



(59) 



(60) 



On the other hand, Fi(z) can be related to the electric field E(r) generated by the electrons 
at the edge of the disk: 



dV 

2poFi{z) = --^ \r^R= E{r)r^R 



(61) 



The divergency of E(r) at r=R is removed when one considers the existence of the edge 
widthness 

EiR) = 2po In (7^) , (62) 

where 7 is a constant that depends on the precise way the electronic density goes to zero at 
R. Thus, 

L 



= 271 p, 



L{L-1) 



1 , / i?\ , A 1 
- /ra 7— + 1 - > 

2 \a ^^2m-l 



(63) 
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If L>>1 but still N>>L, that expression can be further elaborated since 

E^~i(C+tai)+ta2, (64) 

where C is the Euler constant 0.5772... Introducing the wave vector q = L/R, and keeping 
only the leading q-term we get 

Es = uj{q) = ^j2poq\n{qo/q) , (65) 

where go = and /3 = 0.964. 

Result Eq. ([65| ) can also be obtained using classical hydrodynamics It is worthwhile 
to notice that the above ci;(q) differs from that obtained within the edge-magnetoplasmon 
model (Refs. [|TiHT6[): 




which is ^2/3 times that of the two-dimensional plasma frequency. We would like to stress 
that to get the ^Jq ln(go/g) dispersion relation, it is crucial to take into account widthness 
effects in the electronic density. These effects are important in quantum dots, where the 
number of electrons usually is much less than the number of ions, and in any case the electron 
density has a non negligible edge region. 



VI. NUMERICAL RESULTS 

We have applied the method of Section IV to dots made of N= 6, 20, 30, 42, and 56 
electrons. We have taken the values of g*,e and m* indicated after Eq. (|l|), and N+=125, 
i.e., a positive density of ~ 4 ■ lO^^ions/cm^. This corresponds to a GaAs disk of about 1000 
A radius modelling the positive background. 



The g.s. of the dots has been obtained using the CDFT of Ref ||2T]]. We have checked that 
we reproduce their results when we use the same V+. For large dots and intense magnetic 
fields, sometimes one has to face severe convergence problems in the solution of the Kohn- 
Sham equations. Rather than a deficiency of the numerical algorithm, we consider it as a 
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consequence of a inherent characteristic of the system under study, namely the existence 
of a very dense s.p. energy spectrum. To overcome it on one hand, and to carry out 
the calculations under conditions closer to the experimental ones on the other hand, we 
have found it convenient to approach the description of the collective spectrum as the low- 
temperature limit of the results obtained from the finite temperature generalization of the 
formalism of Section II (see Ref. [^]), and of the Kohn-Sham equations (see for example, 
Ref. [^). Consequently, the numerical results we discuss below have been obtained at a 
temperature T ~ 1-2 K. A comparison with several cases in which the T=0 calculation is 
easy to converge, allows us to state that the small temperature we use does not influence 
the results here presented. Thermal effects on the collective spectrum of quantum dots will 
be described in details elsewhere. 

Figures |l|-H show the electronic densities corresponding to dots with N=6, 30 and 56 
electrons for B=0, and for B=5 tesla, respectively. Figures |]-| display the B-dependent, 
multipolar spectrum of the same dots up to L=4. This is the interesting region where the 
crossing between uoj^l and uj^li branches may occur and has been experimentally observed. 
The energies are drawn in meV, and the magnetic field in tesla. These figures show that 
as the uj-L energies go to zero, they may reach a value comparable to the electron-hole s.p. 
energy difference at a rather moderate B-value. When this happens, the collective state 
lies within the s.p. excitations region and experiences a strong Landau damping, loosing its 
collectivity and eventually being washed out. This is one of the reasons why the uj+l branch 
has been experimentally observed up to higher B- values than the uj_i branch [^]. 

The structures appearing along the L>1 curves are due to drastic changes in the particle, 
kinetic and paramagnetic current densities of the dot arising from the effect of B through 
the exchage- correlation energy, which has a profound influence on them (see Fig. 3 of Ref. 



21|). The structures in the uj±i branches roughly correspond to values of B at which the 
total g.s. spin has a minimum. For example, for N=6 we have that 2S is equal to 2 at B=l 
T, to zero at B=2 T, and to 2 at B~3 T. For N=56, one has that 2S=13 at B=4 T, 7 at 
B~5 T, and 8 at B=6 T. For the N=6 dot, the rising of the ooj^l^i curves at B~5 T is due to 
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the full alignment of the electron spins. No similar risings show up for N=30 and 56 because 
for them, the alignment occurs at B-values higher than those displayed in the figures. 

It can be seen from Figs. |^-|^ that the crossing between cu+i and uj^l branches does 
not follow a clear N-systematics. We have also plotted in Fig. ^ the E3(Ql) energies for 
L=l and 4 (dashed lines). The scaling energy is reproducing uj^l to within 10-20 %, since 
the tU-L branch is contributing to m3(Qi) with the same weight as the oj+l one, and the 
negative B-dispersion energies are going to zero rather slowly. 

We also display in Fig. || (dashed lines), the collective energies obtained for a parabolic 
potential whose c^o has been fixed to 5.6 meV in order to reproduce the uji energy at B=0. 
It is worth to notice that this value does not equal the one which fits the Coulomb potential 
generated by the 1000 A radius disk charged with N_|_=125 ions near the origin, which is 4.4 
meV, and consequently, ujq has to be interpreted here as an effective parameter to reproduce 
the dipole energy at zero magnetic field. For L=l, to the scale of the figure both calculations 
coincide. It may be seen that the B-slopes of the uj±l branches for L=2 to 4 are roughly the 
same for the parabolic and disk confining potentials. 

At B=0, Figs. |^-|^ show for N=6 and 56 respectively, the different relative contributions 
to m3 coming from kinetic and Coulomb energies, as a function of L. The Coulomb energy is 
also decomposed into e-e and dot-electron (+e) energies. These figures show that for small 
L-values, m3 is dominated by the (+e) contribution, the kinetic and (e-e) ones being of 
minor importance. However, for a fixed N the kinetic contribution is eventually taking over 
the Coulomb energy contribution. It occurs at an angular momentum Lcr which increases 
when N increases. Since Q^, behaves as ~ r^, the higher the L, the more external its 
influence on the electronic density. It means that for large enough L, it just acts on the 
s.p. wave functions having the bigger 1-angular momentum and radial quantum number n. 
Consequently, it no longer generates collective but s.p. excitations. As collectivity has its 
origin in interparticle effects, it has been argued in Ref. [|13| that Lcr roughly represents 
the largest L collective mode the system can sustain, since for L > the restoring force 
represented by ms is basically determined by an independent particle property like the 
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kinetic energy. This criterion yields Lcr ~ 4, 6 and 8 for N=6, 30 and 56, respectively. 

We represent in Fig. || the energy of modes with L=l to 4 at B=0, as a function of the 
number of electrons in the dot. The ratios oj^/uji and UJ2/UJ1 have, for N>20, average values 
1.68 and 1.37 respectively, instead of 3 and 2 as corresponds to the harmonic oscillator 
sequence. It is apparent from that figure the N-dependence of ojl. As we have already 
discussed, only for the dipole mode in the case of a parabolic confining potential one has 
uoi = ujq. The deviation from that rule becomes more and more important as the number 
of electrons in the dot increases. This has been already discussed in Ref. |T^. The N- 



dependence of ul^i is more complex, and it is due to the interplay between kinetic and 
Coulomb energy contributions to the excitation energy. It can be easily understood from 
Eq. (pHD, which shows that the contribution of the kinetic energy term is important for small 
N and increases with increasing L. For fixed L, the Coulomb energy eventually dominates, 
and ujl>i depends on N in a way similar to the dipole mode. 

Finally, we have also studied the collective spectrum of the N=6 dot in the parabolic 
potential with 0^0=5.6 meV, and have found that as expected, the uj^^i energies at B=0 
depend on N. For example, UJ2 decreases from 8.2 to 7.5 meV, tus decreases from 10.9 to 9.0 
meV, and decreases from 14.3 to 10.4 meV when one goes from N=6 to N=56. 

So far, we have presented a systematic study carried out under well defined conditions 
which could render it difficult the comparison with the results of a given experiment, since 
there are several variables that have to be fixed at the experimental values to permit a 
sensible comparison. To end this section, we want to compare some results obtained within 
our formalism with the experimental data of Ref. [Q. This is possible only in part because 
in that work, the L=l and 2 well defined branches have been obtained for a dot made of a 
large number of electrons, N=210. At present, this is too large a value for us to deal with 
microscopically. Let us remind that the s.p. wave functions behave near the origin as r''' 
(see for example Ref. [0), and since the s.p. levels are nondegenerate when Bt^O, large N's 
imply large s.p. angular momenta. Consequently, we have only attempted to describe the 
N=25, R=1000 A dot. It can be inferred from the value of the dipole energy Demel et al. 
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find at B=0, that N+ ~ 28. For this dot, the ratio UJ2/UJ1 they obtain is around 2, a value we 
are unable to reproduce, whereas for the N=210 dot it is ~1.5, in better agreement with our 
systematics. Nevertheless, it can be seen from Fig. ^that cj+i and E3(Qi) nicely reproduce 
the positive B-dispersion branch, which is the only one for which a detailed comparison is 
possible. 



VII. SUMMARY 

In this work we have used a variational approach similar to the one proposed by Feyn- 
mann in the case of liquid ^He, to describe the multipole spectrum of quantum dots. One of 
the merits of the method is that quantum and finite size effects can be taken into account. 
It may be easily applied to dots hosting several tens of electrons. Rather than in the method 
itself, this limitation has its origin in technical difficulties inherent to current microscopic 
approaches to handle a large number of electrons in intense magnetic fields. 

We have presented a systematic description of edge modes up to L=4 in the region of 
interest to describe level crossing at non zero magnetic fields. We have given explicit formulas 
for the mi and m.^ sum rules corresponding to the general multipole operator Eq. (^91) , which 
reduce to very simple expressions in the dipole case, Eqs. ( |A23|) . These sum rules are exact 
within CDFT, and may be of interest to check the accuracy of any detailed calculation of the 



L-mode strength, in a similar way as they are currently used within TDLDA ||3^. Besides 
this practical application, it is worth to notice that there are few studies in the literature 
of an ma sum rule corresponding to a physical situation where time-reversal invariance is 



violated (see for example Refs. p^ , |35| for the three dimensional polarized electron gas). 



For large neutral dots at zero magnetic field, we have shown that the classical hydrody- 



namic dispersion law for edge waves co'(q) ~ In (qo/q) holds when quantum and finite size 
effects are taken into account. Finally, we have also shown that in the case of a parabolic 
potential, the dipole mode can be exactly solved, yielding the well known classical formula 
for u±i. The exactness of the dipole collective spectrum was stressed in Refs. Here, 
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we have obtained it in a different way, and have gone a step further exphciting the structure 
of the I ±1) collective states. This has allowed us to get a variational solution valid for any 
axially symmetric lateral confining potential. 
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APPENDIX: 



In this Appendix we give some hints about how to derive results Eqs. (^^-^), and the 
sum rules mi and corresponding to the external one-body operator Ql of Eq. (p9D, which 
are given in Eqs. (^) and (0). The details of the method we use here can be found in 
Refs. 0,0,0,0. 

We fix our attention on the operator O^^ of Eq. (0); the operator Ol^ is handled 
similarly. The double conmutator in the numerator of Eq. (|^ can be decomposed in three 
pieces: 

[0+L, [H, = al {[Qt, [H, Ql]] + ih+L [H, Qi] - [Qt, [H, Ri]) 

+ hlL'[Rt[H,RL]\} (Al) 

Splitting the Hamiltonian Eq. (0) into a one body term Hq and a two-body term V=X]j<j 
V (I Ti — Tj I), and using that Ql is a local operator which commutates with V one gets: 

[H, Ql] = [Ho, Ql] = -iLRl + \luj,Ql , (A2) 

where we have used the conmutation relations 

[T, Ql] = -iLRl , [L„ Ql] = LQl . (A3) 

From Eq. ( [A^ ) and 

[Qt,P] = 2tLQU , [L^,P]=P , (A4) 
it is then inmediate to derive for the first two terms of Eq. ([Al|) the following results: 

(0 I [Qt, [H, Ql]] I 0) = 2L^ J r^"^-^ p{r) df = 2m^ (A5) 

(0 I [Rt, [H, Ql]] - [Ql, [H, Rl]] [ 0) = -2i{uo,m^ + ALT) , (A6) 

where 
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r = 1(0 I Rl] I 0) = (L - 1) I r2^-S(r) . (A7) 
To evaluate [i?^, [if, i?^]] we first calculate 

[H,Ri] = [H,,Ri] + [V,Rl] (A8) 

The first term gives 

N . ■ N 

[Ho, Rl] = -i{L - 1) ^ + -Lu.Rl + -cu^'Ql + ^ Y.^^. + ^a3,)K,(ri)gi~i , (A9) 

i=i ^ ^ i=i 

and it is then straightforward to recover the results of Eqs. ( ^OD and (^) for 

t2 n 

m,{T) = -(0 I [Rt,-^{L - 1)^5^-2^'] I 0) (AlO) 
^ i=i 



and 



^2 N 

rhsi+e) = —(0 | , + z9,)K,(r,)gi-i] | 0) . (All) 

i=l 



2 

The remaining two terms in Eq. (|A9|) give: 



2 

which yields the Eq. (^7|) term, and 



(0 I [Rl , ^cu.Rl] I 0) = 2La;,r , (A12) 



{0\[Rt,'-uj'M\0) = l^m, (A13) 



which yields the second term in Eq. (pBI) . 

It remains the problem of evaluating the (0 | [R^, [V,Rl]] \ 0) term. We have done it 
within CDFT by scaling the CDFT g.s. as 



r^) = e^^^ I 0) (A14) 



and then calculating the derivative: 

1 



2 9^2(^1 ^l^)l.=o > (A15) 



where 
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Wf!0Mfil (A16) 

2, J ri — r2 



and 



N 



Pr,{r) = {v I \v) = P + VPl + +V^P2 + ■■■ (A17) 



i=l 



with 



pi = -Lr^-V'(^) (A18) 

P2 = -L\2L - I)r2^-V(r) + -LV^^-y . (A19) 
2 2 



From Eqs. ( |AT^ - |AT9D one gets the result Eq. 0: 



m3(ee) = y(0 | , [V, i?^.]] | 0) . (A20) 

It is seen that the exchange correlation energy does not give any explicit contribution to 
ifi^{ee). However, it affects g.s. magnitudes lilce p(r) and the other densities, thus implicitly 
inffuencing all these quantities. 

The conmutator [O^i^, O^l] of the denominator of Eq. is easily evaluated to be 
(0 I [0+^, 0+l] I 0) = a^(0 I [Ql - tb+LRL , Qt + ih+LRt] \ 0) = 4a^(6+mi + fe^L^E) 

(A21) 

The cubic energy weighted sum rule Eq. ([TOD for the external operator is given by 

ms = ^(0 I [[H, -iLRt + ]^Luj,Qll -iLRl + \luj,Ql] \ 0) , (A22) 

where we have employed Eq. ( [A2D . Using the previous result it is easy to recover the 
expression Eq. (0), which at B=0 reduces to fh^ Eq. (|39D with 7(r)=0. 

Let us finally indicate that for the dipole operator, the mi and ms sum rules for any 
value of B and an axially symmetric confining potential V+ have the simple expressions: 

N 

miiD) = - 

(A23) 
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FIGURES 

FIG. 1. Electronic densities in (ag)"^ as a function of r in ag, for dots with N=6, 30 and 56 
electrons and B=0. 

FIG. 2. Same as Fig. 1 for B= 5 tesla. 
FIG. 3. U!±L energies in meV as a function of B in tesla for L=l to 4, corresponding to N=6. 

FIG. 4. Same as Fig. 3 for N=30. The dashed lines are the E^^Ql) energies for L=l and 4. 

FIG. 5. Same as Fig. 3 for N=56. The dashed lines represent uj±l>i obtained using a parabolic 
potential instead of the one generated by the disk. 

FIG. 6. Decomposition of m3(QL) into kinetic (solid line), total Coulomb (dashed line), 
(e-e)-component (dash-dotted line), and (-l-e)-component (dotted line) contributions as a function 
of L for the N=6 dot at B=0. The lines are drawn to guide the eye. 

FIG. 7. Same as Fig. 6 for N=56. 

FIG. 8. Zero magnetic field ujl energies in meV for L=l to 4, as a function of N. 

FIG. 9. io±i (solid lines) and E3(Qi) (dashed line) energies in meV as a function of B in tesla, 
for a dot of R=1000 A, N+=28, and N=25. The values of e, g* and m* are those of GaAs given in 
the text. The points are experimental results taken from Ref. [6]. 
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